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I INTRODUCTION 


1.1  The  jnethoda ^and  results  of  reference  (l^are  applied  to  the  calcula- 

tion of  the  relative  energy  distribution  among  the  frequency  component?  of  a 
shock  wave  at  various  distances  Aram  the  explosion  of  a one-half  pound  charge 
In  sea  water.  The  mathematical  details  are  sinoarixed  in  Appendix  I.  pjjuali- 
tatively,  the  method  used  and  approximations  involved  may  be  described\ae 
follows i \ 

(a)  At  distances  relatively  close  to  the  charge,  the  shock  wave  is 
represented  as  an  exponentially  decaying  wave  with  a mathe- 
matically discontinuous  front. 

(b)  As  the  shock  wave  travels  outward,  viscous  effects  cause  the 
selective  attenuation  of  higher  frequency  components  in  the 
manner  described  empirically  by  Liebermann  (2).  This  produces 
a change  in  the  shape  of  the  shock  wave,  rounding  off  the  front 
and  decreasing  the  relative  energy  content  in  the  higher 
frequencies. 

(c)  Owing  to  finite  amplitude  effects,  the  profile  of  the  shock 

wave  is  continually  spreading,  resulting  in  a continually  in- 
creasing time  constant.  Also  the  rounding  of  the  front  by 
viscous  attenuation  is  partly  opposed  by  the  effect  of  finite 
amplitude  in  restoring  higher  frequency  components  and 
sharpening  the  front.  These  effects  are  treated  by  means  of 
various  approximations  described  in  reference  (1).  The 
essence  of  the  method  lies  in  the  calculation  of  an  "equivalent 
viscous  distance",  % , which  represents  the  distance  that  the 

shock  wave  would  travel  in  the  Absence  of  finite  mnplitude 
effects  in  order  to  be  subjected  to  the  same  effective  attenu- 
ation of  higher  frequency  components  as  it  actually  experiences 
in  travelling  the  distance  ft  under  the  influence  of  the 
superposed  effects  of  viscosity  and  finite  amplitude. 


II  FREQUENCY  SPECTRUM  OF  THE  ATTENUATED  WAVE 
AT  VARIOUS  DISTANCES  FROM  A 1/2  LB.  CHARGE 

2.1  Figure  1 shows  the  energy  distribution  function  (equation  (17), 
Appendix  I)  plotted  against  frequency.  Energy  flux  is  given  as  db  above 
10-6  ergs/cm?  for  various  distances  from  the  explosion.  The  data  are 
given  in  the  frequency  range  for  which  the  method  of  analysis  is  reasonably 
accurate.  At  frequencies  lower  than  about  0.5  to  1 kc,  equations  (12)  and 
(17)  of  Appendix  I would  be  increasingly  in  error  since  it  is  obvious  that 
the  explosion  wave  is  not  one  of  infinite  duration  and  that  therefore  the 


f 


low  frequency  epectrun  as  given  by  the  exponential  wave  of  equation  (1) 
has  no  physical  weaning.  The  actual  distribution  function  must  curve  down 
and  reach  zero  at  sero  frequency.  At  frequencies  higher  than  30  he,  the 
nethod  of  reference  (1)  introduces  too  severe  an  attenuation  by  neglecting 
the  relaxation  effect  described  by  Liebereann(2) . 

2.2  Figure  2 shows  the  seaa  data  as  Figure  1,  bub  with  the  energy 
distribution  function  plotted  against  distance  for  various  frequencies • 

2.3  Tables  of  the  maearical  results  on  which  Figures  1 and  2 are  based 
are  given  in  Appendix  II. 


APHSHDIX  I 


Mathematical  Development 
1.  The  Uhattenuated  Shock  Wave. 


The  initial,  undistorted  shock  nave  produced  by  the  underwater 
explosion  may  be  represented  as  an  exponentially  decaying  pulse  with  a 
mathematical  discontinuity  at  timet  t ■ o 


O 


9 t >o 
a t<0 


where  A is  the  reciprocal  of  the  tlam  constant  0 and  is  the  peak 
pressure.  The  experimentally  observed  time  constants  and  peak  pressure  i 
be  represented  empirically  > 

a.  >/\ . St  (mt*) 


A/S 


where  # • charge  weight  in  lb. 

R ■ radial  distance  from  explosion  in  ft. 
0 • time  constant  in  sec. 

^ - peak  pressure  in  lb. /in. 2 


The  exponential  shock  wave  of  equation  (1)  may  be  represented  by 
a Fourier  integral! 


p(t)  m 


i 


-u- 


Contolnlng  (1)  and  (5),  it  nay  be  shewn  that 

?e~)  . Z/c*r)*C*  tu>) 

and  that  the  abeolute  magnitude  of  ^ 1st 

|^|  - T. /rg^a^)* 

2,  Energy  fits  In  the  unattenuated  ease. 

The  total  energy  flux  In  the  shock  ease  Is  giren  byt 


whore  pc  Is  the  acoustic  laps  dance  of  the  eater. 
Using  Rancherel's  theorem^ 


and  conbinlng  equations  (7)  and  (9)t 


Evaluating  the  Integral  In  equation  (10)  t 


Em  %*/jlAfC 


The  distribution  function  of  energy  flux  with  respect  to  angular  frequency  1st 

f»y.  l _ 7Z  L 


ec~) 


e?J TfC 


The  normalised  distribution  function  is  thus  given  by: 


£lc->)  - 


■tr  /»v« 


The  phoreieal  significance  of  equation  (33)  is 


presents  the  fraction  of  the  total  energy  flux  contributed 


by  angular  frequencies  in  the  range  between  Af  and  Similarly 

represents  the  fraction  of  the  total  energy  flux  contributed 

< , , 

by  frequencies  in  the  range  between  T,  end  £ • 


3*  Kffect  of  viscous  attenuation. 

Following  the  method  of  reference  (1)  p.  38-59  and  p.  7U-76,  we 
apply  Liebemann's  (2)  frequency-attenuation  results  for  sea  water  to  the 
explosion  pressure  wave  and  introduce  first  order  corrections  for  finite 
amplitude  effects.  The  "equivalent  viscous  distance" , 9C  , is  given  as  a 
function  of  £ in  Pig.  10  of  reference  (1).  Using  Liebemann's  data,  the 
attenuation  parameter,  €t  , is  given  in  terms  of  X i 


a*  7.  //  (/o~*)-?c- 


It  ia  than  sham  that  tha  total  ana  ray  flux,  £ , 
in  the  attenuated  wave  is  given  byt 


where  s\a-J  ia  the  conventional  probability  integral  defined  and 

tabulated  In  various  Tables  of  Functions  such  as  reference  (3). 

The  energy  distribution  function  is  them 


frrpc 

and  the  nonaalimed  energy  distribution  functiont 


Qjr  the  earn  argument  used  above;  j 
energy  flux  in  the  frequency  range 


Tables  of  Energy  Distribution  Functions 


w • Charge  Weight  ■ lb. 

R • Radial  Distance  from  Charge,  ft. 

£%)•  Energy  Flux  Distri^lon  Junction 

4t0*K  S'  . — erga/c*2 

* Ib/in2 

Values  of  4t*and  4*are  given  with  each  table. 


TABIZ  1 


I?  - 1000  ft.j  Jl*«  20.5(106)sec”^}  a?  • 790O(1O“16)|  S - 6.88  Ib/ia* 


- 8 - 
TABLE  2 


R - 2500  ft.  j 13.3<  106)*  - 3U,800(10-16)|  6-  2.1*1*  Ib/i a? 


£%0x 

£%h\  *bov» 

30-6  ergs/o*2 

0.5 

1290 

31.1 

1.0 

563 

27.5 

2.0 

173 

22.1* 

5.0 

29.6 

U».7 

10.0 

7.32 

8.6 

15.0 

3.13 

5.0 

20.0 

1.69 

2.3 

25.0 

1.01 

0.00 

30.0 

0.6$ 

-1.9 

TABU  3 

R - $000  ft.j  A*.  10.K1O6);  a2  - 96,000(10"16)}  1J-  1.12  lb/*»2 


1*01*) x IflJ 

db  above 

erga/oa2 

10^  erga/oBt2 

0.5 

313 

25.0 

1.0 

126 

21.0 

2.0 

36.9 

15.7 

5.0 

6.11* 

7.9 

10.0 

1.1*6 

1.6 

15.0 

0.568 

-2.3 

20.0 

0.291 

-5.1* 

25.0 

0.157 

-8.0 

30.0 

0.088 

-10.6 

7500  ft.;  A*- 


T7ET 

o.$ 

1*0 

2.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 


TABU  1* 

8.35(M^)|  ■ 167,000(10"i6)|  { - 0.70  lb/la? 


•rta/or 


336 

51.7 

3J*.8 

2.1*1 

0.9*8 

0.206 

0.092 

0.01*38 

0.0211 


21.3 

17.1 

11.7 

3.8 

-2.6 


- 10  - 
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